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Abstract 

We consider a perturbed ordinary differential equation where the perturbation is only signifi¬ 
cant when a one-dimensional null recurrent diffusion is close to zero. We investigate the first 
order correction to the unperturbed system and prove a central limit theorem type result, i.e., 
that the normalized deviation process converges in law in the space of continuous functions 
to a limit process which we identify. We show that this limit process has a component which 
only moves when the limit of the null-recurrent fast motion equals zero. The set of these times 
forms a zero-measure Cantor set and therefore the limiting process cannot be described by a 
standard SDE. We characterize this process by its infinitesimal generator (with appropriate 
boundary conditions) and we also characterize the process as the weak solution of an SDE 
that depends on the local time of the fast motion process. We also investigate the long time 
behavior of such a system when the unperturbed motion is trivial. In this case, we show that 
the long-time limit is constant on a set of full Lebesgue measure with probability 1, but it has 
nontrivial drift and diffusion components that move only when the fast motion equals zero. 

1 Introduction 

In systems of multidimensional diffusion processes where both slow and fast timescales are 
present, the limiting behavior of the slow component is an interesting mathematical problem 
with applications in physics, biology and other areas. One possible setting is a system of 
diffusion processes (A e (f), Y £ (t )) that depends on a parameter e representing the ratio of the 
two time scales. More precisely, suppose that X e changes faster and faster in time as e 4- 0 
while Y £ changes on the same time scale for all values of e. Then X £ and Y £ are called the 
fast and the slow component respectively. 

Classical stochastic averaging theory, first invented by Khasminskii ([16]) expanding the 
ideas of Bogolyubov for ordinary differential equations ( 0 ), states that if the fast component 
X s has a finite invariant measure, then as e —>• 0, the distribution of the slow component 
Y £ approaches a limit where the X dependence in the drift and diffusion coefficients of Y is 
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averaged out with respect to this invariant measure. This result was later extended and refined 
by a vast number of authors (see e.g amongst many others). 

In this paper, we are going to study a system where the fast component is a null-recurrent 
process and therefore has no invariant probability measure while the slow motion is a per¬ 
turbation of an ordinary differential equation. Namely, let IT be a fc-dimensional Brownian 
motion and consider the system 


dX £ (t) = e 1 ip(X £ (t),Y £ (t))dW(t), X £ (0) = xo/e, (1.1) 

dY £ {t) = b{X £ (t ), Y £ (t))dt + cr(X £ (t), Y £ {t))dW{t), Y £ { 0) = y 0 , (1.2) 

where all coefficients are Lipschitz continuous, p is bounded away from zero and infinity, and 
the process X £ ( Y £ ) is one (d) dimensional. Intuitively, the fast component X £ (t ) is of order 
0{e~ l ) and therefore the limiting dynamics of the system, as e 0, is governed by the behavior 
of the coefficients for large absolute value of the x-variable. 

In 0~7]> it was shown that if b, <p, and a have certain limits as |x| f 00 in the Cesaro 
sense then (eX £ (t), Y £ (t )) converges in law to a process where the x-dependence in dra¬ 
ma) is replaced by coefficients associated to these limits. More precisely, after the change of 
coordinates (X £ (t), Y £ (t )) = ( eX £ (t ), Y £ {t)), the system ( X £ {t ), Y £ (t )) can be written as 

X £ (t) \ = f 0 \ / (p(e~ 1 X e (t), Y £ (t)) 

Y £ (t) ) V b{e~ l X £ (t),Y £ (t)) ) C a(e~ 1 X £ (t), Y £ (t)) 

It might be tempting to assume that the limit will be given by taking the limit of the coefficients 
for large positive and negative values of x. However, this is not the complete picture as one has 
to include the effects of null recurrent averaging. To see this, note that the diffusion matrix 
can be written as 



^ dW(t) 


A(x, y) 


1 
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W(x,y )\ 2 cpa T (x,y) 
&<p T (x, y) cra T (x,y) 
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i yv T (x,y) \ 
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Let p(x,y ) = l/\(p(x,y )\ 2 and note that if we freeze the y dynamics, this is proportional to the 
sigma finite invariant density of the X £ process. For a function K(x, y), define 


K±(x, y) 


i r u 

X{x>0} • lim - / K(v, y)dv + X{x<o} 
1 1 u-H-oo u Jo 


i r u 

lim — / K(v,y)dv, 

u^-oo U J o 


whenever the Cesaro-limits exist. Using this, introduce the averaged diffusion matrix A{j(x, y) = 
( pAij)±(x,y)/p±(x,y ) for i, j = 1,.., 1 + d and the averaged drift b(x,y ) = (0, {pb)±/p±). 

Theorem 1.1 (Khasminskii, Krylov [T7]). Suppose that <p,cr,b are Lipschitz continuous and, 
for each x, their first and second order derivatives in y are bounded continuous functions of y. 
Assume as well that for some positive constants c\, 02,03 > 0, we have 


d 

ci < \p(x, y )| < c 2 , ^ 2 (A t (x, y))a + bf{x, y)) < c 3 (l + \y\ 2 ) 

i =1 
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for all (x,y) € M 1+ci . Also assume that p,pAij,pbi and all their first and second derivatives 
in y have well defined Cesaro limits as x —>• Too. Then the pair (sX £ (t),Y £ (t)) converges in 
distribution in C ([0, oo), R 1+d ) to the solution of 

d (X°(t), Y°(t )) = b(X°(t), Y\t))dt + yjA(X0(t),Y°(f))dW(t), 
whenever weak uniqueness for this equation holds. 

In general, however, the Cesaro limits in question can be quite different for large positive 
and negative values of x. This implies that the coefficients of this limiting system will be 
discontinuous at x = 0 in the generic situation and the well-posedness of the corresponding 
martingale problem becomes a challenge. Note that the limit being well defined is an assump¬ 
tion in this theorem. For a list of special cases when well-posedness holds, we refer the reader 
to m (see also ED- 

In this paper, we answer the question what happens when the Cesaro limits of a are zero 
while the Cesaro limits of b are equal. Namely, we study the following special case of (ll.l|) - (11.2D . 

dX £ {t) = £- l y(X £ (t), Y £ (t))dW(t ), X £ (0) = x Q fe, (1.3) 

dY £ (t) = [bi(Y £ (t)) + b 2 {X £ (t), Y £ mdt + a{X £ (t\ Y £ (t))dW(t ), E £ (0) = y 0 , (1.4) 

where b 2 and a satisfy a certain decay condition in x which guarantee that their Cesaro limits 
for large \x\ vanish so their contribution to the dynamics for large values of |X e | is negligible. 
In this setting, the result of m implies that Y £ (t) converges uniformly in probability to the 
solution y(t) of the ordinary differential equation (ODE) 

y(t) = b\ (y(tj). y(0) = y 0 • 

In fact, it can be shown that this convergence holds in a stronger sense (see Lemma 13.21) . 
Consequently, Y £ (t) can be viewed as a perturbation of y(t ) and it is of interest to describe 
the behavior of their deviation. This system can model, for example, the effect of spatial 
inhomogenity localized around a one co-dimensional surface (e.g. a crack in a solid) where the 
homogeneous structures away from the interface can possibly be different on different sides of 
this interface. This is the situation, for example, in groundwater hydrology when one studies 
diffusion in composite porous materials, see e.g. [20] and the references therein. In our scaling, 
the dynamics transversal to the interface is much faster than the one parallel to it. 

In [22] . the authors considered the simpler case when X £ is a Brownian motion independent 
of the noise driving the slow motion. Namely, the stochastic differential equation (SDE) 

dY e (t) = [&i(E £ (f)) + bAe-'Wfit), Y £ )]dt + aie-'Wfit), Y £ (t))dW 2 (t), Y £ { 0) = y 0 , 

where W\ and W 2 are independent standard Brownian motions, was considered and it was 
shown that the first order correction to the deterministic limit is of order yje. Moreover, a 
functional central limit theorem was proved for the normalized difference e _1 / 2 (y £ (t) — y{t)) 
with the limit identified as a process that is Gaussian when conditioned on W\. This case was 
significantly simpler than (ll.3l) - (ll.4l) as the difficulty arising from the discontinuity at x = 0 
does not appear and there is no feedback from the Y £ process to X £ . This allowed for a 
proof using only elementary stochastic analysis. In the limit, the perturbation acts only when 
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W\ (t) = 0 (i.e., on a Cantor set of times) and the equation for the limiting solution can be 
expressed in terms of the local time of W\ at zero. 

In the current paper, we prove a similar result for the general system (ESD-dHD using the 
martingale method. We show that if a ^ 0, the pair — y(f))) converges to 

a non-trivial process with singular diffusion on the interface which we describe. On the other 
hand if a = 0, then the pair (eX £ (t), £ _1 (y e (t) — y(t))) converges to a nontrivial process with 
singular drift on the interface. Moreover, we show that in the absence of the unperturbed 
system (b± = 0), the interesting limit behavior happens on longer time scales (of order £~ 2 ). 
The long-time limit has both singular drift and singular diffusion in the sense that there is a 
proper diffusion process on the interface, with coefficients that are averages with respect to a 
sigma-finite invariant measure of the fast motion, that is sampled at the random times when 
the limit of the fast motion is in the origin. Once again, these times form a Cantor set and 
the limiting slow motion is constant for almost all times. 

It is worth pointing out that in some sense, the original result of m and the results in this 
paper parallel the first few steps of the program contained in (DEI, Chapter 7) when the fast 
process is positively recurrent (See Section 7). 

The plan of this paper is as follows. In Section [2] we specify the conditions on the coefficients 
in (USD-dEID, describe the limit processes, and state our convergence results. In Section [3l we 
derive the main technical tools needed for our analysis. In Section HJ we show that the families 
of processes under consideration are tight while in Section 0 we prove that all subsequential 
limits converge to the solution of the associated martingale problem. In Section [6j we prove 
that these martingale problems are well-posed and characterize their unique solutions as the 
weak solutions of certain SDE-s. In Section [71 we mention some open problems. 
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2 Assumptions and main results 

For convenience, we will make the change of coordinates (X £ (t),Y £ (t)) = (eX £ (t),Y £ (t)) and 
consider the system 

dX £ (t ) = ip(e~ 1 X £ (t), Y £ (t))dW(t ), X £ (0) = x 0 , (2-1) 

dY £ (t) = [6i {Y £ (t)) + b 2 (s~ 1 X £ (t), Y £ (t))]dt + a(s- 1 X £ (t), Y £ (t))dW(t ), T £ (0) = y 0 , (2.2) 

where W(t) is a fc-dimensional Wiener process, tp(x, y) is a A;-dimensional vector field, b\(y) and 
b- 2 (x, y) are d-dimensional vector fields, and c t(x, y) is a d x fc-matrix. In this way, X £ (t) is a 
one-dimensional process and Y £ (t) is a d-dimensional process. 

We make the following assumptions: 

1. ip, & 2 , and o are Lipschitz continuous and b±(y) is a twice continuously differentiable 
d-dimensional vector field with bounded derivatives. 
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2. Decay conditions hold in x: 


d 

S(x) := sup \b 2 (x,y)\ R d G D 1 (M), cf 2 (x) = supTr aa T (x,y) = sup V' cr, 2 (x, y) € L^R). 
2/GR 2/6R 2/GR ^j = 1 

(2.3) 

3. JW is non-degenerate. More precisely, if \ip(x,y)\ 2 = Yli=i ^fbDy), we require that there 
are numbers 0 < c\ < c 2 < oo such that 


ci < | <p(x, y )| 2 < c 2 V(x, y) G R 1+d 


4. \/\ip(x, y)\ 2 has Cesaro-limits in x, i.e. the limits 


lim — 

x —»±oo x 



1 

W{u,y)\ 2 


du 


exist. Define 


a±{x, y) = X{x>()} ( lim - [ 
L - 1 yu— kx> u J o 


0 l ^( D 2/)| 2 


dv 


-l 


+ X{x<0} Jm 


' u Jo 


-°°uj 0 \ip(v,y)\ 2 


dv 


-l 


(2.4) 

Occasionally, we will use the notation a + (y) to denote the restriction of a± to {(x, y) € 
R d+1 : x > 0} (similarly we will use a_(y)). 

It follows from Assumption 1 that the ODE 


— = bi(y(t)) y( 0) = 2 / 0 , 
dt 


which serves as the unperturbed part of the slow motion, has a unique globally defined solution. 

Before stating our main result, let us introduce the construction which we use to describe 
the limit processes. Let H + = {(x, y) : x > 0} and = {(x, y) : x < 0} be two half-spaces in 
M 1+d and let C± be second order differential operators acting on functions over H± respectively. 
Denote the interface by Hq = {(x,y) : x = 0}. Also let /3i(y), otij(y) € C(M d ), i,j = 1, ...,d, be 
such that ( oiij(y )) is uniformly elliptic and such that the SDE 


dJJ(t) = j3(U(t))dt + y/ a(U(t))d,W(t) 

has a unique weak solution. Let p±(y) > 0 with p + (y)+p_(y) = 1. Let Co(M 1+d ) be the Banach 
space of bounded, continuous functions / : M 1+c( —> M decaying at infinity endowed with the 
supremum norm. 

We construct an operator C on C'o(M 1+rf ) as follows. Let its domain V(C) C C'o(M 1+d ) 
consist of / G C'o(M 1+d ) such that 

• The restriction of / to H + , H and along Hq are smooth. 

• d y J'(x. y) G C' 1 (M 1+<i ) for i = 1,..., d. 

• £+/(0+, y) = C-f(0-,y), for all y G R d . 
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• d x f(x,y) has right and left limits as x —> 0 and the gluing-boundary condition 

d rl 



(2.5) 


is satisfied for all (x, y) € R 1+d . 


For any / € D(£), let Cf(x,y ) = C±f(x,y ) for x € H± and define Cf(0,y) by continuity. 

This is the adaptation of a special case of the well-known Wentzell type boundary con¬ 
ditions (see [3D], [3TJ or the survey [29|1. In general, it is hard to prove that the closure of 
such an operator is the generator of a continuous strong Markov process Z(t ) = ( X(t),Y(t )). 
Nevertheless, we will show that when p+(y ) = p~{y ) = 1/2 and C± has a special form, C does 
indeed generate a Markov semigroup on Co(R 1+a! ). In fact, in the cases of interest we construct 
the corresponding Markov process probabilistically. However, we believe that the closure of 
such an operator C is the generator of a Markov semigroup in greater generality. 

Intuitively, the process Z(t) can be described as follows. Inside the half-spaces H + and 
17_, it coincides with the diffusion generated by C± up to the point when it hits Hq. On 
Hq, Z{t ) follows a diffusion process for an infinitesimal amount of time, after which it leaves 
Hq in the direction of H± with probability p±. The set of times at which Z(t) € Hq has 
Lebesgue measure zero and consequently the displacement due to the effect of the boundary is 
singular in time (the support of the increments is a Cantor set in R). To make this intuition 
rigorous and obtain a characterization of Z(t), let us recall the definition of the local time of a 
semimartingale (see Theorem 7.1 in EH). 

Definition 2.1. The symmetric local time of a semimartingale S with quadratic variation (S) 
is the unique nonnegative random field 


L = {L s (t, x); (f, x) € [0, oo) x R} 


such that the following hold: 

1. The mapping (t,x) —» L s (t,x) is measurable and L s (t,x ) is adapted. 

2. For each x € R, the mapping t —>• L s (t, x) is non-decreasing and constant on each open 


interval where S(t) x. 


3. For every Borel measurable f : R —> [0, oo), we have 



a.s. 


R 


4- L s (t,x ) is a.s. jointly continuous in t and x for x f 0, while the one sided limits exist 
at x = 0 and 
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Note that instead of (12.61) . we could have choosen an arbitrary convex combination. Other 
common choices include L s (t, 0) = L s (t, 0+) (right local time) or L s (t, 0) = L s (t, 0—) (left 
local time). For our current purposes, however, the symmetric local time satisfying (12.61) is the 
most suitable as it makes our formulas more transparent. 

Assume that Z{t ) is the unique solution of the martingale problem given C for appropriate 
j3, a and p±. It is natural to conjecture (see also e.g. mm) that it is also the unique weak 
solution of an SDE involving singular terms on the interface x = 0. Namely, let 

C± = ^A(x, y)d xx + B(x, y)d y (2.7) 

be the diffusion operator on i7_ U H + , where d y is the d-dimensional gradient with respect to 
the y variable and A(x, y ), B(x, y) are matrix-valued functions that are uniformly elliptic and 
sufficiently regular on both halfspaces H 1 . Assume A(x,y) has a jump discontinuity at x = 0 
(B(x, y) may or may not be continuous at x = 0). Then Z(t) is a weak-solution to the following 
SDE 

dZi(t) = y/A(Z\(t), Z 2 (t))dW(t) + [p + (Z 2 (t )) - p-(Z 2 (t))\dL Zl (t), Zi(0) = x 0 , (2.8) 

dZ 2 (t) = B{Zi(t), Z 2 (t))dt + (3(Z- 2 (t))dL Zl (t) + yJa(Z 2 (t))dV Zl (t ), Z 2 ( 0) = y 0 , (2.9) 

where y7 is the matrix square root, L Zl (t ) = L Zl (t, 0) is the local time of Z\ at 0, IE is a one¬ 
dimensional Brownian motion. V Zl (t) is a continuous martingale whose quadratic variation 
is (V Zl ,V Zl )t = $ijL Zl (t) and the cross-variations are ( V Zl 1 W)t = 0 for all i = 1 ,...,d. It 
is natural to define the weak solution to such an SDE as a quadruple (Z 2 , Z 2 , W, V Zl ) such 
that IE is a one-dimensional Brownian motion and V Zl is a continuous martingale satisfying 
the above variation relations and that (Zi,Z 2 ) solves (I2.8I) - (I2.9D in the integral sense with 
(IE, V Zl ) taken as the driving noise. In this way, Z\ diffuses like a skew-Brownian motion (see 
pTj ) modulated by Z 2 (t), and Z 2 (t) — B(Zi(s), Z 2 (s))ds moves only when Z\(t) = 0, because 
L Zl (t) (and therefore V Zl (t)) is constant when Z\ ^ 0. 

Because of the discontinuous coefficients and the coupled nature of the above equations, 
the well-posedness of such an SDE is in general unclear. We will show, just as one would 
expect based on the classical situation, that the weak solutions of dSBII-lQ are in one to 
one correspondence with the solutions of the martingale problem associated to the operator L 
as discussed above under fairly general circumstances. However, we do not need this general 
result for the proofs of our main results and therefore defer the proof to the appendix. 

We are ready to state our main result on the asymptotic behavior of the process (X e , Y £ ). 
As we mentioned before, Y £ (t) converges to y(t). Our first main result concerns the rate of 
this convergence. Let C([0, oo), M fc ) be the space of all continuous M fc -valued functions on the 
halfline endowed with the metric 


f yi 1 BUPte[ 0 ,n]l/(*)~g(*)l 

P ,9 “{ 2n 1 + su PtG[0,n] 1/(0 - 5(0| 

It is well known that weak convergence of probability measures over this space is equivalent to 
weak convergence of measures in CdO, T],M fc ) with the supremum norm for all T > 0. Let us 
introduce the normalized deviation process C e (0 = £^ 1 ^ 2 (Y £ (t) — y(t )). 
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Theorem 2.2. The triple {X £ {t),y{t)X £ {t)) converges in law in C([0, oo), M 1+2d ) to a Markov 
process which is the unique solution of the martingale problem associated to the operator C 
described above with 

1 Q Q 

C ± = 2 a±(x ’ y) d^ + b i(y) • fyj + (dyhix, y)- Q- — 

where a± is as in \2-4\ ) and 

° ^ 

(w) (x ’ y)dx ) ' 

Moreover, there is a one dimensional Brownian motion W and a continuous martingale V x ° 
such that (V x ° ,V x °)t = 5ijL x ° and (Vi,W) = 0 for each i = 1 and such that the limit 

process of the pair ( X £ ,( £ ) satisfies the inhomogeneous equations 

dX°(t) = \/a±(X°(t), y(t))dW(t ), (2.10) 

d(°(t) = d y bi(X°(t), y(t)) • C °(t)dt + (. x ,y(f))dx^dV x °(t), ( 2 . 11 ) 

with (X°(0), C°(0)) = (cco, 0 ) where yf- in \2.11\) denotes the matrix square root, and d y bi(y ) is 
the derivative tensor of the vector field b\ at y £ W l , i.e ( d y bi(y))ij = d(bif / dyT 

Remark 2.3. Note that once the processes X° and y are constructed, the solution of \2.11\) 
is explicitly given by the variation of parameters formula 

C°«) = J‘ (^ 1 ) <x, MXtz) avx\s). 

As is apparent from this theorem, the drift part of the perturbation (b 2 ) does not contribute 
to the deviations of order yje. On the other hand, the following result suggests that 62 does 
play a role in fluctuations of order e. 

Theorem 2.4. Let a = 0 and ( £ (t) = e~ l {Y £ {t) — y(t)). The triple (X £ (t), y(t), ( £ {t)) converges 
in law in C([ 0 , 00 ), M 1+2d ) to the unique solution of the martingale problem associated to the 
operator L described above with 

1 [) Q 

c ± = 2 a±(x ’ v) qt2 + h i(y) • + (W®. v) ■ 0 • ^ 

and 

p±{y) = 1 / 2 , m = ( 0 , r 

\ J —OO 

Moreover, there is a Brownian motion W such that the limit process of the pair (X £ ,£ £ ) solves 
the inhomogeneous equation 

dX°{t ) = yj a±(A°(s), y(s))dW{t) A°(0) = x 0 , 

d(°(t) = d y bi{X°{t), y{t)) • C \t)dt + (^J (^ 2 ) y(t))dx^ dL x \t) C°(0) = 0 , 

where in \2.11\) denotes the matrix square root, and d y b\(y ) is the derivative tensor of the 
vector field b\ at y € i.e ( d y b\{y))ij = d(b\ f /dy ] . 


b 2 

M 2 


(x, y)dx 


a(y) = 0. 


P±(y ) = 1/2, /%) = 0, a(y) = 

















Remark 2.5. In this case, the variation of parameters formula implies 


C(t) = 


f 

Jo 


o f*(dyb 1 (X 0 (r),y(r))dr 


f 


2 ^ (x, y(s))dx ) dL x (s). 


/o \J -oo \M 2 

When bi(y ) = 0, we have y{t) = 0 and it follows from Theorem 12.21 that sup te[0T ] |W(f)| 
converges to zero for every fixed T in probability. On long time scales (of order e -2 ), however, 
the process converges to a limit with both singular drift and singular diffusion on the interface. 
Let X s (t ) = eX £ (e~ 2 t) and Y £ (t) = Y £ (e~ 2 t). This pair satisfies the SDE 

dX £ (t) = (p(£- 2 X e (t), Y £ (t))dW(t), X £ (0) = ex 0l (2.12) 

dY £ (t ) = e~ 2 b2(E~ 2 X £ (t), Y £ (tf)dt + £~ 1 a(e~ 2 X £ (t), Y £ (t))dW(t), T e (0) = yo, (2.13) 

with a Brownian motion W different from the one in (12. 1 D - (I2.2D . 

Theorem 2.6. The pair ( X £ , Y £ ) converges in law in C([0, oo); M 1+rf ) to the unique solution of 
the martingale problem associated to the operator C with 

r 1 i 

C±= 2 a±{X ' V) d?- 


and 


p±(y ) = 1/2, /%) = 


f 


(x, y)dx, a(y ) = 


f 


GO 


\t\ 2 J ’ 7-oo VM 2 

Moreover, the limit process is the unique weak solution of the equation 


(x, y)dx. 


dX\t ) = Ja±(X°(t),Y<\t))dW(t), 


(2.14) 


dY°(t) = 


b 2 

l^l 2 


( x, Y°(t))dx 


dL x ° (t) + 



(T(T ± 




(x, Y°{t))dx 


dV x \t ), 


(2.15) 


with (X°(0), y(0)) = (0, 2 / 0 ) where L x °, W, and V x ° 


are as in (USD and m - 

As one can see from (12.151) . the singular drift and diffusion on the interface are both 
significant on these longer time scales. Also, a careful look at the explicit formula for j3{y) and 
a(y ) reveals that they are essentially averages with respect to a u-finite invariant measure of 
the fixed y process dX y {t) = ip{X v (t),y)dW(t, ) (with density c/\ip(-,y )\ 2 for any c > 0). In this 
sense, Theorem 12.61 is an averaging result where the singular term on the interface follows the 
averaged dynamics of the slow component with respect to the invariant measure of the fast 
motion. Due to the null recurrence of the fast dynamics, however, the natural timescale for 
this averaged process is not proportional to t (as it is due to the ergodic behavior of fast-slow 
systems in the positive recurrent case). Instead, it is proportional to the amount of time that 
the fast component asymptotically spends in the essential support of b and a, or in other words 
to the local time at x = 0. Since our fast motion is roughly Brownian, this time scale is of 
order \ft. 

Note that even though the theorems above are stated for a flat interface, they are applicable 
to situations when the interface is a more general hypersurface. We also remark that we chose 
the initial conditions xq, yo somewhat arbitrarily. The modifications one needs to make to 
include the case when xo,yo can depend on e differently are straightforward. 
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3 Auxiliary lemmas 

First we derive an estimate on certain integral functionals of the fast process m 
Lemma 3.1. For any p > 1, there exists C p > 0 such that for all if € ZdfR), all stopping 

v/2 

times 0 < ti < T 2 with Erf < oo, and all e > 0, 

< C p \if\ P Llm E\T 2 -n|2. (3.1) 


i r 

^ J n 


l X £ (r))dr 


Proof. If / € C* 2 (R), then by Ito formula, 

f(8- l X £ (T 2 )) - Re^X^i ri )) = (3.2) 

= ~J f'(£~ 1 X £ (r))tp(£~ 1 X £ (r),Y £ (r))dW(r) + ^ J f"(E~ 1 X £ (r))\ip(£~ 1 X £ (r),Y £ (r))\ 2 dr 

almost surely. In fact, by the Meyer-Tanaka formula, the above expression is valid even if f" 
is not continuous, but is only L 1 . We want to construct / such that fix) = \if(x)\. Since 
if € L^R), the function 

T(.x) = / \if{y)\dy 

J —OO 

is well defined, continuous and bounded (0 < ’P(x) < |V , Il 1 (R))- Then we can define 

fix) = f 'P iy)dy. 

■Jo 

This function is Lipschitz continuous because 


\f(x 2 ) -fix i)| < 


rx 2 

J Xl 


V(y)dy 


< I^IlHK)I X 2 - Xi\. 


(3.3) 


Now, because ci < \tp(x, y)\ 2 , 


| ifie L X e (r))\dr < ce 


i r 

2 


e X £ ir))\\ipi£ X £ ir),Y £ ir))\ dr, 


where here and in what follows, c is some (possibly different from line to line) positive constant 
independent of e. Plugging this into (13.21) . 

1 f T2 

- / \ifi£~ 1 X £ ir))\dr < 

^ J n 

< c (eifiE - 1 X £ ir 2 )) - fi£~ l X e in))) - [ ne- l X £ ir))vi£- l X £ ir),Y £ ir))dWir)^ 
By the elementary inequality \a + b\ p < 2 p ~ 1 (|a| p + \b\ p ), 


+ E 


if is 1 X £ ir))dr 


< 


c(eeP | fi£~ l X £ iT 2 )) - fi£~ 1 X £ iT 1 ))\ I 




f 

J T-\ 


\P(e 1 X £ ir))(pis 1 X £ ir),Y £ ir))dWir) 


_ —1 V£ 
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By & 

£\f(£- 1 x e ('my) - fie-'x^in))] < 

< £ I'0 |l 1 (M)I £_1 ^ £ ( t 2 ) - £ _1 X £ (ti)| = \lf\ L l(R) 


f 

j T1 


ip(e 1 X £ (r),Y £ (r))dW(r) 


and therefore, 

1 f T2 


i n 
^ J T\ 


i/)(£ 1 X £ (r))dr 




+ cE 


r T 2 

J n 


</>(e 1 X e (?’), Y £ {r))dW{r) 


+ 


^(s- 1 X £ (r))ip(£~ i X £ (r),Y £ (r))dW{r) 


.-1 V£ 


Each of these expectations can be bounded by the BDG inequality, and we see that 


i r 2 


i/}(£ l X £ {r))dr 


'n 


— C l'^lz, 1 (R) 


r r 2 


+ c 


e j 

r*T2 


| ip(s 1 X £ {r),Y £ {r))\ A dr 


p/2 


+ 


'h 2 (e~ 1 X £ (r)) | ^(e~ 1 X £ (r), T £ (r)) | 1 "dr 


— 1 V£/ 


f n 


p/2 


All of these integrands are bounded because |y>| < C2 and 'h(x’) < |Vi.Li(R)- Therefore, we can 
conclude that 

1 /*T9 p 

-jl^ ^ , £ 

izJ(R) 1 " "" 


i r 
- ./1 


i/>(£ 1 X £ (r))dr 


< CM P Ll(m E\r 2 - n\ E 2. 


as required. 


□ 


Next we show that the slow motion (12.21) converges to the unperturbed system in a strong 
sense. 

Lemma 3.2. For every p > 1, there exists a constant C p such that 

sup E\Y £ (t) - y(t)\ p < C p max{e p T p / 4 ,e p / 2 T p / 2 }e^ TP_1 . 

te[0,T] 

Proof. Note that 

Y e (t) ~ V(t ) = [ (bi(Y £ (s)) - bi(y(s))ds + [ b 2 (e _1 A £ (s), T £ (s)) ds+ 

Jo Jo 


a 


{£~ i x £ (s), y £ ( S )) dw(s) = i\(t)+m)+my 


By the Lipschitz continuity of b\ and Jensen’s inequality, 


E\im p < T p ~ L Lipihf / E|y £ (s) - y(s)\ p . 

Jo 

On the other hand, recall that b{x) = sup ye]R d \b 2 (x,y)\ and note 


sup E|/f(t)| p < E 
te[0,T] 


\P 

b(e~ 1 X £ (s))ds 1 <C p £ p T p / 2 , 
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where in the last inequality we used Lemma 13.II with t\ = 0 and 72 = T. Finally, it is easy to 
see that the scalar quadratic variation of I 3 is 

(h)t= I Tracr T (s~ 1 X e (s),Y £ (s))ds< f a 2 (e~ l X £ (s))ds, 

Jo Jo 

(where we recall cr 2 (x ) = sup ygR d TV era 1 ( x , y)) and therefore by the Burkholder-Davis-Gundy 
inequality and Lemma 13. II we have 


/ r T \ p/2 

sup E|/f(f)| p < C P E[ a 2 (e~ 1 X £ (s))ds ) <C p e p/2 T p/4 . (3.4) 

te[0,T] \J 0 J 


The result now follows from Gronwall’s lemma. 


□ 


Next, we recall Lemma 3.5 from CZI, which explains why the Cesaro limits of describe 
the limiting behavior of X s . 


Lemma 3.3. For any f G C(M+ x M x R d ) and any T < 00 , we have 

f f(s , TG(s), T £ (s)) - a±(X £ (s), Y £ (s)))ds 

Jo 


sup 

t€[0,T] 


0 , 


(3.5) 


as e 4- 0 in probability where a± is defined in (El Moreover, \3.5\) holds for every family of 
stopping times {r £ } e < eo in place of T such that {TG(t £ )} £ < £0 is tight. 

Before proving Lemma 13.31 we establish the following fact about Cesaro limits. 

Lemma 3.4. Assume that g : M —>• R is bounded and has the Cesaro limit 

1 


g := lim 

x —>-+oo x 


g(Od4- 


Then for any f G L 4 oc (R), and R > 0, 


lim sup 

£— ^0 ire[0,-R] 


* r - / 

Jo Jo 


g I f(Qd£- I f(Og(£ 1 0df 


= 0 . 


(3.6) 


Proof. First, we consider the case where / = 1 (see m). In this case, for i£l, 

j (g - g(£- l O)dfi = x ^g — ~ J giOd^j =■ xa . 


(3.7) 


Because g is the Cesaro limit of g, for any fixed x, the above expression converges to zero as 
e —>• 0. This convergence is uniform in |x| < R because for |x| < R, 


xa 


(-) 

— 

[ (g - g(.£ 1 0 )dfi 

V £/ 


Jo 


< min 


2x\g\ L oo,R 



(3.8) 
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Next, we consider the case where / £ C 1 . In this case, by integrating by parts, 



/(0(5-0(e 1 0 )<% 


f(x)xa - j ( ~ £ j d f. 


(3.9) 


Therefore, 


sup 

0<x<R 


f 

JO 


f(0(g-g(s 1 0 )<% 


< sup 

xa 

[ 7 ) 

0<x<R 




i\f\cH0,K\) Q- + R) (3.10) 


and it follows from (13.811 that the convergence is uniform. 

Finally, for a general / £ ^i' oc (®0 and R > 0, we can approximate / in T 1 ([0, R]) by a 
sequence {/„} C C 1 ([0, R]). Then for n € N, 


sup 

o <x<R 



f(0(g - g(£ l O)df 


< 2\g\L™([-R,R])\f-fn\L 1 ([-R,Rl)+ SU P 

0 <x<R 



(0(g - g(e 1 0)d£ 


By first choosing n large and then e small we can make the above expression arbitrarily 
small. □ 


Proof of Lemma, 1 3. ‘A Very similar to the proof of Lemma 3.5 in [IT] . we only include a sketch 
of the proof. For e > 0, define the auxiliary function u £ (t , x, y) to solve 


\ip(s l x , y)\ 2 u e xx (t, x, y) = fit , x, y)(\ip(e 1 x, y) | 2 - a±(x, y)). 
u £ (t, 0, y) = u%{t , 0, y) = 0, lim £ ^ 0 sup| x |< R \u £ (t, x, y)\ = 0. 


(3.11) 


By dividing by |<^>| 2 and integrating, we see that for x > 0 (the situation for x < 0 is 
similar), 

f x r x fit £ v) 

u £ x (t,x,y) = j f(t,£,y)d£- a + (y) J y )\2 d ^ 

By Lemma 13.41 this converges to zero as e —> 0 uniformly in \x\ < R if and only if a+iy) is 
the reciprocal of the Cesaro limit of . By integrating in x once more, we show that u £ also 
converges to zero as e —>• 0 uniformly in |x| < R. The proof concludes by applying Ito formula 
to u £ (t, X £ {t), Y £ (t)), identifying the integral in (13.51) as the Ito-correction term, a standard 
localization argument and estimating the resulting other terms. The last statement of the 
lemma follows similarly using that X £ (t £ ) is in a large enough compact set for all e < £o with 
high probability. □ 


Finally, we state the version of the Ito-Tanaka-Meyer formula that is relevant to our situ¬ 
ation (a special case of [23]). 

Lemma 3.5. Let (Zi, Z 2 , W, V Zl ) be a weak solution of <f 2.S\) - (2.g\) with p + = p- = 1/2 and 
assume that f £ (7(M 1+rf ) has continuous first and second order derivatives in y and for fixed 
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y € M d , /(., y) G C' 1 (M_ t _ U R_). Then 


f(Z\(t), Z 2 (t)) = f(Zm, Z 2 m + f d^fiZ^s), Z 2 (s)) y /A(Z 1 (s), Z 2 {s))dW{s) 

J o 

^ ^ rt 

+ yz d y J(Zi(s), Z 2 (s))yJa(Z 2 (s))dV z \s ) 

i=i Jo 


+ l 

d i d 

+ X] ®yif( z i( s ')’ z 2<,s))(3i(Z 2 ) + - ^2 dy iVj f(Zi(t), Z 2 (t))aij(Z 2 (t)) 


i =1 


P+OT>))/ X (0+, Z 2 (t)) - p_(Y(s))f x (0-,Z 2 (t)) 

d 

2 ' 


*J=1 


dL z \t ) 


1 f l 

+ 7T / dxxKZxis), Z 2 (s))A(Zi(s), Z 2 (s))ds + 

^ io 



d yi f(Zi(s ), Z 2 0))-B(Zi0), Z 2 (s))ds, 


where used the symmetrized x-derivative 

d sym /{x y) = + 

Note that the symmetrized derivative was necessary by our choice to use the symmetric 
local time. 


4 Tightness 

In this section we are going to show that the laws of the families in Theorem 12.21 Theorem 12.41 
and Theorem 12.61 are tight in the space C([0, oo); R 1+rf ). To prove tightness, we will show that 
there exists c>0, p>0, q>l such that for all T > 0 and 0 < s < t < T, and e > 0, 


E\X £ (t) - X e (s)| p < c\t-s\ q , 


E|C £ (t) - C £ (s)r < C\t - s| 9 

and in the long-time case 

E\X £ (t) - X £ (s)\ p < c\t-s\ q , 
E\Y £ (t)-Y £ (s)\ p <c\t-s\ q . 
All of these results are consequences of Lemma 13.11 


Theorem 4.1. Let (X £ ,(^ £ ) satisfy the conditions of Theorem 12.21 or Theorem \2.f\ Then 
the family of laws of the pairs {(X £ , C £ )}e>o ore tight. Let ( X £ ,Y £ ) satisfy the conditions of 
Theorem \2J^ then the family of laws of the pairs {(X £ , T e )} e> o are tight. 

Proof. We prove tightness in the context of Theorem l2.21 as the tightness for the other systems 
can be proved analogously. In this case, 


X £ (t)-X £ (s)= I ip(£- 1 X £ (r),Y £ (r))dW(r) 
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By the BDG inequality and the fact that \tp\ 2 < C2, 


E| A' e (t) - X £ (s)| p < c 2 \t - s\ p/2 . 

As for (f £ , because Y £ (t) = y(t ) + y/e( £ (t), 

Cm - CM = [ + dr + e - 1/2 J\^ x%r)Mr) + 

+ £ -1/2 J (7(s~ 1 X e (r),y(r) + cCC(ry)dW(r) = If(s, t) + /f(s, t) + /f(s,i). 
We estimate /f(s,f) and I^(s,t) by using Lemma (3T1 By (12.31) and Lemma PTT1 


E|/f(s,t)| p <E 


r 1/2 / b(e~ l X £ (r))dr 


< ce 1 / 2 ^ — s| p / 2 , 


(4.1) 


and similarly by the BDG inequality, 


//•t \ p/ 2 

(e _i y ( j 2 (e“ 1 X £ (r))drJ < c|t - s| p/4 . (4.2) 


E|if(s, t)| p < cE 
Finally, if p > 1, Lemma 13.21 implies 

E\Il(s,t)\ p < ^11611101^1)1* “»I P 


(4.3) 


We choose p>4, q = p/A and tightness follows from (14.ID . (14.21) . and (14.31) and the Kolmogorov 
criterion above. □ 


As a result of tightness, for every subsequence e n —>• 0 there exists a e nk such that 
( X £n k , £ En fc) converges in distribution. In the next sections we characterize the limit and 
show that it is the same process for all such subsequences. Consequently, the entire sequence 
converges in law to that unique limit. 


5 Convergence to the martingale problem 
5.1 Freidlin-Wentzell type result 

We are going to show that the processes in Theorems 12.21 12.41 and 12.61 each converge to a 
solution of a martingale problem. In all three cases, the main tool in proving convergence 
is a Freidlin-Wentzell type result (®. see also [12]). Assume we have a family of (1 + k)- 
dimensional processes Z £ . Depending on the situation, Z £ will be either Z £ = (X £ ,y,( s ) 
(k = 2d) or Z £ = (X £ ,Y £ ) (k = d). We are going to use the notation E,(-) to the expectation 
when Z £ has initial point z = ( x , y, Q or z = (x, y) respectively. Recall that H± = {z : ±x > 0} 

Theorem 5.1. Let Ci be second order differential operators on Hi, i = {+,—} and let Di C 
C£°(Hi) be a set of test functions that have bounded derivatives of all order. Let t £ = inf{t > 
0 : \X £ (t)\ < 1(e)} for some 1(e) | 0 as e | 0. Assume that there is a Ao > 0 such that for each 
A > Aq, the following assumptions hold. 
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1. We have for every f E Di 


T S 


r e 

f(Z e (r £ )) - f(z) + / e~ xt (A f(Z £ (t)) - C t f(Z £ (t))) dt 

Jo 


= 0(k(e)) (5.1) 


for some k(e ) | 0 as £ | 0 uniformly for z E Hi. 

2. For i E {+,—}, there is a Ui t \ E Di such that 

£iUi,\ = A u i>x 'iz € Hi, |x| < 1 

and such that u± i \(Q, y) = 1 and ±d x u± t \(0, y) < c < 0. 

3. There is a 5(e) f 0 as e j. 0 such that 5(e)/1(e) —>• oo, 5(e)/k(e) —>• oo and 


sup E z e xt x { - S (e),S(e))(X £ (t))dt ->• 0 
ZgRi+d Jo 


(5.2) 


as e | 0. 
f. Define 

9 5 = inf {t > 0 : |X £ (t)| > 5(e)}. (5.3) 

There are p+(y), p~(y) > 0 such that p + + p~ = 1 and 

P z(Z"(0 s ) E Hi) —$■ pi(y) 
as e | 0 uniformly in z E M 1+rf , |.x| < 1(e). 

5. There are /3i(y),cxij(y) fori,j = 1 such that 

\w.(Yf(0 & ) ~ Vi) ->• A;(y) (Y £ (0 5 ) - yi)(Yf(9 5 ) - y 3 ) -»■ c^(y), 

as e | 0, uniformly for |x| < 1(e). Moreover, we have supm < ^ e - ) E\Z £ ( 9 s ) — z| 3 = o(5). 

Then for every f E Co(M 1+d ) such that f\ui E Di and the boundary-gluing conditions (USD 
hold, we have 


ess sup 


E 2 



[a f(Z £ (t)) - Cf(Z £ (t ))] dt - e- xt °f(ZDM)\F t0 


—> 0, 


as e | 0 for every to > 0, A > Ao, where (Ft 0 = cr(Z £ (s); s E [0, to!)- In particular, for all such 
functions, every weak limit of Z £ as e | 0 satisfies the martingale problem for L defined by 
£f(z) = £if(z) for z E Hi. 


Proof. By the Markov property of (Z £ ,P z ), it is sufficient to show that 

POO 

A(e, z) = E z e~ xt [X f(Z £ (t)) - Cf(Z £ (t ))] dt - f(z) 0 

Jo 

as e f 0 uniformly with respect to z E R 1+d . 
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Introduce the sequence of stopping times do = 0, (p n = inf{f > a n : |X £ (t)| = 1(e)} and 
cr n _|_i = inf{f > 4> n : |X £ (i)| > <5(e)}. Using this sequence, we can decompose the above ex¬ 
pression to terms corresponding to downcrossings (on [a n , </>„]) and upcrossings (on [c/) n , (J n+ i\) 
respectively. More precisely, 


00 r r<t>n 

A(e, z) = J2 E z e~ x ^f(Z £ (c/ n )) - e~ XrT " f(Z £ (a n )) + / e ~ xt [X f(Z £ (t)) - Cf(Z £ (t))\ dt 

—n L J CTn 


n =0 
00 


+ 


^ fO-n+l 

+ e ~ XtT n+i f (Z £ (a n+i)) - e~ x ^f(Z £ (</> n )) + / e" At [A f(Z £ (t)) - Cf(Z £ (t))\ dt 

n= 0 J<t>n 

Note that by the strong Markov property, we can bound the downcrossing terms by 


Eft 

71=0 


e "E Z = Z £ (p n ) 


-A T e 


f(Z £ (r £ )) - f(z ) + [ T e~ xt [X f(Z £ (t)) - C±f(Z £ (t))] dt 

Jo 

OO 

<0(k(e))J2 Ee 


< 


A CTn 


71=0 


(5.4) 


where we used (EH). To bound the sum on the right hand side, first note that 

OO OO 

J2 Eze~ Xan < 1 + ^ E z e- Xa ^E Z e (<Jn _ l)e - XTe (5.5) 

71=0 71=1 

Let |x| — $(£) and apply ED to Ui t a given by Condition 2. to get 

E z e~ XT " u itX (Z £ (r £ )) - u^ x (z) = 0(k(e)). 


This implies that 

E,e- AT£ = E 2 e- Ar£ (l - u itX (Z £ (r £ ))) + u itX (z ) + 0(k(e)) (5.6) 

Since |Af £ (r £ )| = 1(e), we have |1 — Ui tX (Z £ (r e )\ = 0(l(e)) whereas the condition on the derivative 
of Ui t a normal to the interface implies Ui )X (z) < 1 — C5(e) for some C > 0 since \z\ = 5(e). 
Plugging these back to (15.61) gives 

sup E z e _Ar£ < 1 — C'5(e) 

\z\=8{e) 

for small enough e for some C > 0. This, together with (|5.5D . implies that E z e~ XrTn = 

0(1/5(e)) which in turn implies that the bound in (|5.4I) converges to zero. 

To bound the uperossing terms, first note that the sum of the integrals f^ n+1 can be 
bounded by 

poo 

11 (A - C)f | |ooE z / e- xt X (-8(E)As))(Z £ (t))dt -> 0 

Jo 

as £ | 0 uniformly in 2 ; by (15.21) . where we used that / has bounded derivatives of all order. 
The remaining terms can be written as 

OO OO 

E 2 (e- A ^ +1 - e x ^)f(Z £ (a n+1 )) + ]T e~ x ^(f(Z £ (a n+1 )) ~ f(Z £ (cj> n ))). (5.7) 

71=0 71=0 
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where the first term can immediately be bounded by 


I I/I loo sup E z / Ae xt x ( _s(e),6(e))(Z e (t))dt^0. 
zeR 1+d J 0 

For the second term in (|5.7|> . we can use the strong Markov property to write 

OO 

Y, E * e_A ^ E *'=ZH4>n) ( RZ £ {0 S )) - f(z )) . (5.8) 

n =0 


Using \x'\ = 1(e) and Taylor’s formula, we get 


E z , (f(Z £ (9 6 )) - f(z)) = E z ,[f(Z £ (9 s )) - my')] + 0(l(e)) = 

= 5P z (Z £ (9 5 ) € H + )d x f(0+,y') - 6P z (Z £ (9 s ) € H_)d x f(0-,y')+ 

d i d 1 

+ 6 Y, d if (°> y')^EAY £ (9 5 ) - y[) + 5 Y W(0, y')-^ E AYf(9 s ) - y' i )(Y £ (9 5 ) 

i=l i,j =1 

+ o(5), 


2/j)+ 


where we used the last statement of condition 5. Now using condition 4., the rest of condition 
5. and the gluing conditions, we get that this entire expression is o(5) uniformly in \x'\ = 1(e). 
Noting that 

OO OO 

YE z e~ Mn < Y e ~ Mn = ^( 1 /^), 

71=0 71=0 

115.81) converges to zero and the proof is finished. That every weak limit of Z £ converges to 
the solution of the martingale problem of C follows exactly as in [9] after noticing that if the 
Laplace transform of a function is identically zero for sufficiently large A then the function 
itself is identically zero. □ 


We are going to need some estimates on the exit time 9 s . 

Lemma 5.2. 

E z 9 5 < C5 2 , E z (9 5 ) 2 < C5 4 . 

Moreover, there is a 5 q > 0 such that 6 £ (0, <5q] and |a;| < <5 implies 


cos (Cu5) 

Proof. Let W be a Brownian motion such that X £ (-) and x + W(J Q |<^| 2 (s)ds) have the same 

distribution and therefore so do 9 ^ (the hitting time of 5 by W) and J ® |<£>(.s)| 2 d,s. The results 
follow from this, |<,c| 2 € [ci,C2] and well-known Brownian formulas. □ 
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5.2 Martingale Convergence for the long-time asymptotics 

Now we show that the system of equations in Theorem 12.61 satisfies the conditions of Theorem 
15.11 with Di = As b\ = 0, we denote b = 62 in this case. Let Z £ = (X £ ,Y £ ) solve 

I dX £ (t ) = v(e~ 2 X £ (t), Y £ (t))dW(t), 

\ dY £ {t) = £-' 2 b(£~ 2 X £ (t ), Y £ (t))dt + £-V(£- 2 A £ (f), Y £ (t))dW(t). ' ’ * 


Lemma 5.3. Lei / S L; * = -,+. If we choose 1(e) so that 1(e) —> 0 and £ 2 l{£) —> + 00 , then 
(15.11) is satisfied for Z £ = (X £ , Y £ ) and Cf(x,y ) = \a±{x,y)f xx {x,y). 

Proof. By Ito’s formula and a standard localization argument, 


AT(e, z) = E 2 


e~ Xre f(Z £ (r £ )) - f(z) - f 

Jo 




xnz‘m - <# 


= £" 2 E 


+ -£- 2 E 


f Q e_A * ( XJ bi(e~ 2 X £ (t), Y £ {t))f yi (Z £ (t))J dt 

rr £ ( d \ 

J 0 e ~ M I J2 (^ T fifi£- 2 X £ (tf Y £ (t))f yi yfiZ £ (t))\ dt 

+e _1 E z J e_At J P T )i(£~ 2 X £ (t),Y £ (t))f xyi {Z £ {t)) S j dt 

+l-E z [ T e“ At (|^| 2 ( £~ 2 X £ (t),Y £ (t)) - a±(Z £ (t))) f xx (Z £ (t))dt, 
z Jo 

:=/f+Jf + /f+ Jf. 

We estimate each term separately. Note that 

I{<c£~ 2 E z [ e~ M b{£~ 2 X £ {t))\\f\\ c idt, 

Jo 

then because t £ = inf{f > 0 : |X e | < i(e)}, we have 


I\<ce WfWc^z / e b{£ X £ (t))l {{ x Ht)l>m} dt 

J 0 

00 rk -\-1 

<c£~ 2 \\f\\ c i^2e- Xk E z / b(£~ 2 X £ (t))l {l x e(t) \ >m} dt. 

k =0 Jk 

By Lemma 13. 11 there exists a constant c, independent of k and e such that 

rk+l ^ /■ 

£' 2 E 2 / 6(£" 2 X e (i))l { |_Y S(t) | > ; (£)} (ii < c / S(x)dx. 

{|x|>e _ 2 Z(e)} 


Because 6 € L^M) and £~ 2 1{e) — >• + 00 , it follows that If — >• 0 uniformly in z G Hi. The 
analysis of If is analogous to that of If and we can conclude that If —>• 0 unifomly in z G If,;. 
To estimate If, we recall that ||^||^ < C 2 , and therefore a little manipulation shows that 

\Ii I < e _1 c||/||c 2 E 2 [ e~ xt d(£~ 2 X £ (t))dt. 

Jo 
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By Holder’s inequality for the integral E z fj e xt ■ dt, 



Therefore, If —>• 0 unifomly in z £ Hi by the same arguments used for If. We show that 
If —>• 0 by applying Lemma 13.31 □ 


Lemma 5.4. For each A > 0, and i £ {+, —}, there is a £ Di such that 

7 }ai(z)d xx u i} \(z ) = A u it \(z) 'iz£H i , \x\ < 1 

and such that ^±^( 0 , y) = 1 and ±d x u± ! \(0, y) < c < 0. 

Proof. Set 

~r / ^ 'a; 

«±,a(z, ?/) = e V . 

This function satisfies all the requirements. □ 


Lemma 5.5. Let 5(e) —> 0. Then (15.21) holds for Z £ . 
Proof. First we observe that 


f*k -\-1 


„QO ^ 

/ e~ Xt X(.-6(e),8{e))(X e (t))dt < V e~ Xk E z I X(-S(e),S(s))i X£ (t))dt. 

Jo k=0 Jk 

By Lemma 13.11 there exists c > 0 independent of e and k such that 

pk -\-1 

E 3 / X(.-5(e)Me)){X e (t))dt < c|X(-S( £ ),5( £)) |li < «5(e). 

Jk 

Therefore, for any z € M 1+rf , 

roo 

-Xt 


e A X(- 5 ( £ ), 5 ( £)) (A' £ (t))(it < ^ _ e _ x 5(e) -£ 0 


as e f 0 . 

Lemma 5.6. For the system Z £ , p + = ^ in the sense that 


lim P (Z £ (0 5 ) £ H + ) = lim P (Z £ (6 5 ) £ HJ) = 


1 


£—>-0 £—>-0 2 

uniformly in z £ M 1+rf ; \x\ < 1(e) as long as 1(e)/5(e) —>• 0. 


□ 


(5.10) 


20 





Proof. Let 


P Z (X £ (6 5 ) E H + ) = f £ (x,y) 


Since the exit only depends on the x component, and the effect of ip is a time change that 
does not effect the exit probabilities, it is easy to see that the corresponding Brownian formula 
holds, i.e. 

f £ (x) = + ^( £ )) 

2 o(e) 


and if 1 (e) <C 5(e), 


lim sup 

£_>0 |x|<Z(e) 


p x (X £ (8 s ) € H+) - 


1 


0 , 


as e 4- 0. The part for iL_ follows by complementation. 


□ 


Lemma 5.7. Fori,j = 1 let 

r°o u. 

U l 


Pfy) = 


(era 1 \ 


(x, y)dx, Qiijiy) = f V T.i 9 IJ ( x,y)dx. 

J — < 


J —oo |(£)|2 V ^ J _ oo 1^2 

Then if 5(e), 1(e) —> 0 such that 5(e)/1(e) —»• oo, and 6 s = inf{t > 0 : |X e (i)| > 5}, 


1 


1 


-E Z (Y[(9°) - Vi ) -> f3fy), -E z (Yf(6°) - yi )(Yf(0 o ) - yf) -> c^(y) 


as e —> 0 /or |x| < 1(e). Moreover, 


sup E|F/(e 5 ) - y\ 3 = o(6). 
M <l(e) 


Proof. Observe that 


1 


1 


-E z (Yf(9 d ) - Vi ) = ~-^E z J ^ bfe~ 1 X £ (s), Y £ (s))ds 
By Lemma 13.11 and the Lipschitz continuity of b, and the fact that 


(5.11) 


I h(e 2 X £ (t ), Y £ (t)) - bfe X(t), y)\ < 2 b(e X(t))x{\ e -ixe (t) \> K } + c\Y £ (t) - y\X{\xf<e 2 K}, 
we have 


5e 2 


1(5, e) = xyy sup 
|®|<Z(e) 

2 

6 


r e s 


r0* 


E z [ bfe 2 X £ (s), Y £ (s))ds - E z I bfe 2 X £ (s),y)ds 
Jo Jo 


(5.12) 


< ^VE6t\\t {x > K }b\\ L i + ^ E z 1 {|x^(,)|<A^2}l^ e ( s ) - V\ ds - 
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Using E ,9 s < C5 2 for all | x| < Z(e), we see that the first term can be made less than rj/2 for 
any 77 > 0 if K is chosen sufficiently large. Fix such a K and note that by Cauchy-Schwarz 
inequality, the second term in 115.121) can be bounded from above by 



sup | Y £ (t) 
te[ o, 0 6 ] 


[° s 

y I / !{! Xe(s)\<Ke 2 }d s 


< 



E z sup 

V mo, 0 s ] 


Y £ (t) 




^|x e (s)|</-s:£ 2 }^' s 


1/2 


(5.13) 


where the second term is easily seen to be less than cbEO 5 ) 1 / 2 < c6 by Lemma 13.11 On the 
other hand, (a + b ) 2 < 2(a 2 + b 2 ) and the BDG inequality imply that 


sup | Y £ (t) -y | 2 < 

, te[O,0 5 ] / 


(5.14) 


< 2E 


u 


e ~ _ \ 2 
b(e~ 2 X e (s))ds ) + 2cE 


' i r 6 

Y 2 Jo 


a(e 1 X £ (s))ds 1 < cE z 6° < c5 z , 


where Lemma 13.II and Lemma [5.21 were used in the last two inequalities. By (15.131) . and ()5.141) . 
the second term in (15.121) is less than c5 < rj/2 and consequently 1(6, e) < ij for sufficiently 
small 6 . Since rj was arbitrary, (15.111) implies that 


-E z (Yf(9 5 ) - Vi ) = 


1 f° S 


bi(e 2 X £ (s), y)ds + o(l). 


By the properties of the local time, the above expression is equal to 


£ 2 6 Ez J _ 0 o M 2 


bi , _ 


1 


(e 2 u,y)L x \u,6 5 )du = ^-E z 


r k 

Loo m : 


(u, y)L x \e 2 u, 6 s )du 


where L x \u,t ) is the local time of X £ at u. We calculate by the Tanaka formula that for 
fixed x, 


\e z L xe (e 2 u, 0 5 ) = -E ( I X s (0 s ) — e 2 u| — \x — e 2 u 
o o V 


Because X £ (9 5 ) = ±5, the above expression is bounded by 

6 + Ixl 


and if \x\ < 1(e) <S 6 (e), 


lim -E Z L X "(eu, 6 s ) = 1. 
£—>-0 0 


Therefore, it follows from (15.21) and the dominated convergence theorem that 


-E Z (Y £ (0 S ) - yd 


L 


bi_ 

-oo l^l 5 


(x, y)dx. 
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Similar calculations show that 



and a further straightforward application of Lemma 13.11 yields 


sup E|lf (fl 5 ) - y | 3 = 0(5 3 / 2 ) = o(5). 

\x\<l(e) 


□ 


Set 7 > 0, l(e ) = e 2-7 and 6(e) = max{e 2-27 , e _7 sup^^^) |fc(e)|} where k(e) is the con¬ 
vergence rate in Lemma 15.31 Lemmas 15.3115.71 and Theorem 15.11 together with tightness imply 
that every subsequence of (X s , Y £ ) has a further subsequence converging to a solution of the 
martingale problem associated to C. In the next section we will argue that this martingale 


problem is well posed which establishes the convergence of the entire sequence. 

5.3 Convergence to the solution of the martingale problem in the normalized 
deviation case 

In this section, we show that the assumptions of Theorem 15. II are satisfied by Z £ = (X £ ,y,C : £ ) 
with Z e (0) = z = (x, y, Q under the assumptions of Theorem 12.21 or Theorem 12.41 We start 
with a simple growth estimate on ( £ . 

Lemma 5.8. For any stopping time r > 0 with Et < oo, 

(1) under the conditions of Theorem \2.f\ we have 



(5.15) 


(2) under the conditions of Theorem \2J^ we have 



0 <t<T 


(5.16) 


Proof. Recall that under the conditions of Theorem 12.41 we have 


dC(t) 


1 


[h(y(t) + e( £ (t )) - bi(y(t))\ + ^ b 2 (e 1 X £ (t), y(t ) + e( £ (t )) dt, 


£ 


which implies 


X 


L 


■t 


|C £ (0| < l?l + Il^ilie 1 / IC e ( s )|rf s + 


e 1 b(e 1 X £ (s))ds. 


By Gronwall’s lemma, 
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and therefore by the Cauchy-Schwartz inequality, 


E z sup |C £ (t)| 2 < C (E,e 4 H bl llci-) 1/2 f|^| 4 +E Z ( f e^ke^X^dt 
te[0,r] V 7 y \J 0 

and Lemma l3.il and Lemma 15.21 yield (1). 

To prove (2), note that under the conditions of Theorem 12.21 

d( £ (t ) = -^=[bi(y(t) + y/eC, £ (t)) - b x (y(t))] + ~^=b 2 {e~ l X £ (t), y(t ) T \/eC £ <X>) dt+ 
V£ ve 


1 


H 7=°" (e 1 -XT e (t), y(0 + \feC{t)) dW(t) 
yJ e 


and one can easily see as above that 

E z sup |C £ (t)| 2 < C (E z e 4|l?,l|| c ir 
te[0,r] V 


1/2 


\f\ 4 + e 2 E z 


e 1 b(e 1 X £ (t))dt ) + E 


(i V " (£ " 


X £ (s))ds 


1/2 


where we used the BDG inequality in the last term. (2) once again follows by applying Lemma 
13.II twice and Lemma 15.21 □ 


Lemma 5.9. There is a Ao > 0 such that for A > Ao, if we choose 1(e) so that 1(e) —» 0 and 
e~ l l(e) —>• Too, then <EU) is satisfied with some rate k(e) for Z £ (t) = (X £ (t), y(t), ( £ (t)) with 

£f(x, y, 0 = )^a±(x, y)fxx(x, y , o T h(y) ■ f y (x, y, 0 T (d y h(y) ■ 0 ■ f^x, y, 0. 

Proof. We only prove the lemma under the hypothesis of Theorem 12.21 the other case being 
almost identical. Using Ito’s formula, 


„-Ar e 


/’ 

J 0 


f(Z £ (r £ )) - f(z) T / e~ xt (\f(Z £ (t)) - C t f(Z £ (t)))dt 


= h(z,e) + I 2 (z,e) 


where Ii(z, e) contains terms that can be dealt with exactly as we did with their counterparts 
in Lemma 15.31 and therefore 

sup \h(z, e)| —>• 0. 

zeHi 


On the other hand, 


h(z, e) = \-E z [ 
1 J o 


e Xt fxx(X £ (t),y(t),( £ (t)) 


b\(y(t) + \fe( £ (t)) - h(y(t)) 


- (d v bi)(y(t))( £ (t) 


dt, 


and because b\(y) E Cd(WL d ), we have that 


roo 

\h(z,e)\<Cy/e e~ xt E z \( £ (t)\ 2 dt, 

Jo 


(5.17) 
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where the constant C depeds on \f\c? and the second derivatives of b\. By (I5.16p . we have for 
all t > 0 

E,|C £ (t )| 2 < Cife C2t 

for some C\ , C 2 >0 and therefore choosing Ao = C 2 makes the integral in (15.171) converge for 
every A > Ao and the result is proved. □ 

Lemma 5.10. For each A > 0, and * = +,—, there is a Ui t \ € Di such that 

b~'i ^7.A A* 

Proof. The operators C± are both generators of non-degenerate diffusion processes Z = ( X , y, Q 
in H±. Dehne 

t = inf{t > 0 : X(t) = 0}, 

then it is well known that 

Ui, X (z) := E,e- AT (5.18) 

has the property we desire. That the requirement on the derivative is also fulfilled once again 
follows from a time change argument and the corresponding Brownian formula. □ 

The statements of the following lemma are proved completely analogously to those in the 
previous section. 

Lemma 5.11. Property (15.21) holds for Z e . Also, p+ = p- = in the sense that 
lim P 2 ( Z e (r e ) € H + ) = lim P 2 (Z e (r E ) £ H-) = ^ 

£->0 e-s-0 2 


uniformly for (x, y, z) € M 1+2d , |x| < 1 (e). 
Lemma 5.12. Fori,j = 1 ,..,n, let 


/ %(y ) 


' I-ooW (X ’ y)dX 
0 


ifcr = 0 , 
otherwise, 


Oiijiy) 



(cra T )ij 

M 2 


(x, y)dx. 


Then if 5(e)/1(e) —>• oo ; 

^E 2 (Cf( 0 5 ) - &) -> Pi(v), ^E z ((i(0 S ) - 6 )( (j(8 5 ) - £i) ->• otijiy) (5.19) 

as e —>• 0 uniformly in \x\ < 1(e). Moreover, 

E Z \(U0 S ) - £\ 3 = o(5). (5.20) 

Proof. We only prove the case a = 0, the other one being similar. First note that by (15.151) 
and Lemma E2 there is a 5$ > 0 and M > 0 such that 

E z sup \( £ (t)\ 2 <M, (5.21) 

te[O,0«] 
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(5.22) 


whenever 5 < do- Now write 
|E,«f(« 5 ) - (,) = 


= k E * 


r 1 r e 

J [h(y(s) + s( £ (u)) - h(y(s))] dt + —E z J b 2 (e X e (s), y(s) + e( £ (s))ds 


The second term can be dealt with exactly as in (15.121) and (15.131) and it can be shown to 
converge to 

'° 0 b 2 


m = 

uniformly in |a:| < 1(e) once one notices that 


/ 


-oo M 


; (x,y)dx 


E z sup \y(t) + e( E (t) - y\ 2 < CE Z sup \y(t) - y\ 2 + e 2 E z sup |C £ (t )| 2 = o(l) 
teio.e 5 ] <e[0,6»-5] te[0,6»<S] 

by (15.211) . a standard Gronwall estimate on y(t), and Lemma 15.21 On the other hand, the 
integral in the first term in (15.221) is immediately less than or equal to 

/ \ V 2 


WbiWcE 


\( £ (t)\dt < 


9 5 ■ sup |C £ (t)| 

ie[O,0 5 ] 


< E z sup |C £ (t)| 2 V5~ 2 E z (d s )i 

V te[0,eO / 


which converges to zero by (15.151) and Lemma 15.21 

The second claim in (15.191) and (15.21)1) can be proved similarly. 


□ 


If we set 7 > 0, 1(e) = e 1-7 and <5(e) = maxje 1 " 27 , e _7 fc(e)} then Lemmas 15.9115.121 and 
Theorem [5T] with (y, 0 in place of y show that the possible subsequential limits are all solutions 
of the martingale problem associated to C, and the convergence will be established once we show 
well-posedness in the next section. 


6 Uniqueness and Markov property of the martingale problem and charac¬ 
terization of the solution. 

In this section, we prove that the martingale problems associated to the operators in Theorem 
E 2 Theorem 12.41 and Theorem 12.61 are well-posed and the unique solutions satisfy the stated 
SDEs. 

Our main tool will be the following result which is a special case of a result of Ethier and 
Kurtz which we reformulate for our current purposes. 

Theorem 6.1 ([?( Theorem 4.4.1] Uniqueness of the Martingale Problem). Let C be lin¬ 
ear and dissipative on C'o(M (i+1 ). If D(C) = Co(M d+1 ) and there exists a A > 0 such that 
(A — C)(D(C)) = Co(M d+1 ), then the martingale problem is unique for C. That is, if Z, Z are 
processes in CffO, oo), K 1+d ) and Z( 0) = Z( 0) in distribution and for any f € D, 

f(Z(t)) — f Lf(Z(s))ds and 

Jo 

are martingales with respect to the filtration generated by the coordinate process then Z and Z 
are equal in distribution. Moreover, if there is a solution Z, it is a Markov process corresponding 
to the semigroup on C'o(M d+1 ) generated by the closure of L. 


f(Z(t)) - [ Lf(Z(s))ds, 
Jo 
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In all three cases, we probabilistically construct a solution of the corresponding stochastic 
differential equation and conclude by the Hille-Yosida theorem that the generator of the asso¬ 
ciated Markov semigroup satisfies the hypothesis of Theorem 16.11 This together with Lemma 
13.51 implies that the measure generated by this process is the unique solution of the martingale 
problem associated to C. 


6.1 Long-time case 

By the concrete form of /3 and a (as in Theorem 12.611 . there is a weak solution (Yq. Wf) to the 
equation 

dY 0 (t) = P(Y 0 (t))dt + y/ a(Y 0 (t))dW 2 (t). 

Take an independent Brownian rnoton W\ and consider the time changed process 

Z(t) = (X(t).Y(t)) = (W\(t), Y(j(L Wl (t, ()))) 

where L Wl (t, 0) is the local time of IV\ at zero. Finally, set 

1 


t(s) = 


J 


-du 


'o a±(Wi(u), Y(u)) 
and let s(t) denote its inverse. Let Z{t) = Z(s(t)). 

Lemma 6.2. Z(t) =: (X°(t),Y°(t)) is a weak solution of \2.1f\ )- \2.15\) . 

Proof. First note that if we introduce V Wl (t) = W2(L Wi (t. 0)) then Y satisfies 

Y(t)= f p(Y(s))L w \ds, 0)+ f V 'a{Y{sj)dV w \t) 

Jo Jo 

Note that V Wl is a continuous martingale and 


*(£f 


a i:j (Y(s))d(Y w \V^)(s) <oo=l 


for all t > 0 and therefore the stochastic integral is well defined for all times. Indeed, it 
follows from the construction that (v ^ 1 , V^ Vl \ = 1 (t, 0 )5ij. We also calculate that by (12.311 


^2 aii 


2=1 


Tl , ( ;; \ x,y)dx < c\a 12 


Therefore, 




J2 / (V ™ 1 , V ^ 1 ) {s) < c\a\ 2 L2m L w \t, 0). 


which is finite with probability 1. Now performing the second time change yields 
X°(t) = Wi(s(f)) 


Y°(t)= / f3(Y°(u))L w \s(du), 0)+ / a(Y°(u))V Wl (s(du)) 
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Note that for any Borel measurable / : R —>• [0, oo) 


ps(t) pt 

f(x)L Wl (s(t),x)dx = / f(Wi(u))du= / f(Wi(s(u)))a±(W 1 (s(u)), Y(s(u)))du 

Jo Jo 


I 


= I f(X u (u))a±(X 0 (u),Y°(u))du = I f(x)L*\t,x)dx, 


x°, 


which implies L" x) = L x °(t,x). This also implies that if we let V x °(t ) = V Wl (s(t)), 
then we have 


(V x ,V X )(t) = (V, Wl ,Y m )(s(t)) = 5ijL Wl (s(t), x) = 5ijL x (t,x). 


This means that 


L 


Y u (t)= I p(Y°(u))L x \du, 0)+ / a(Y"(u))dV*\u) 


x° r 


as required. We also have that the martingale 


Wit) : = 


~L 


o ^ / a±(X 0 (u),Y°(u)) 


Wi(s(du)) 


has quadratic variation t and therefore by Levy’s theorem it is a Brownian motion. Conse¬ 
quently, 

X u (t) = J* Ja±(X°(u),Y%u))dW(.u). 


It remains to show that the appropriate cross-variations vanish. Since (Wi, W 2 ) = 0, it is 
easy to see that (W \, V^ V ' ) t = 0. As s'(t) € [ci, c 2 ], (X°, V x °) t = (W\,V^ Vl ) s ( t) = 0. Moreover, 


(W,V x °) t = 


l 


0 y/a±(X°(u),Y°(u)) 


d(X °, V x °) u = 0. 


□ 

Proof of Theorem I '2. hi As a Markov process, Z(t) generates a strongly continuous contraction 
semigroup Tt on Co(M 1+rf ) by Ttf{z) = E z f{Z{t)). From this definition and Lemma 13.51 it is 
easy to see that the closure of C (denoted by L again) is an extension of the generator of T t 
and hence, by a standard semigroup fact, they are actually equal. Thus, by the Hille Yosida 
theorem, T>(£) is dense in Co(K 1+d ) and 1Z( A — £) = Co(M 1+rf ). Therefore by Theorem 16.11 the 
associated martingale problem is well-posed. It is straightforward to show using Lemma 13.51 
that Z(t) is actually a solution of the martingale problem. This, combined with the tightness 
result of Section 4 and the convergence result of Section 5, finishes the proof of Theorem 

□ 
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6.2 Normalized deviation case 

Similarly to the previous case, we are going to directly construct the limiting Markov process 
which solves the martingale problem associated to C. We are only going to carry out the 
construction, the other steps are analogous. 

Lemma 6.3. There is a process X° which solves 

dX°{t ) = \/ a±(X°{t), y(t))dW(t). 

This X° solves the martingale problem associated to C for functions f that do not depend on 
V or (. 

Proof. Such a result is not completely obvious because of the discontinuity of the coefficients. 
Consider instead the system 


dy(s) 


dX{s ) = dW (s) 

h(y(s)) j 

-*- as. 

a±(X{s),y{s)) 


The pair (X, y) is perfectly well-defined because the formula for y is a random ODE that we 
can consider pathwise in the integral sense. The function b\(y)/a±(x, y) is discontinuous in x 
, but uniformly Lipschitz continuous in y. Classical successive approximation arguments show 
that a unique solution exists. We then define the random time change 

f s 1 

t(s) = / --- dr 

Jo a±(X(r),y(r)) 

and denote its functional inverse by s(f). Then we define 

X°(t) = X(s(t)), 
y(t) = y(s(t)) 


which has the properties that we desire. □ 

The following two lemmas are a straightforward consequence of Lemma 13.51 
Lemma 6.4. 


1. Let (X e ,y,( £ ) satisfy the assumptions of Theorem A2.il Let 

CM) = l e/s °y bl( y (r » dr M j (aa T )(x, y(s))dx^j dV x °(s). (6.1) 

where V x ° = Wo(L x °(t, 0)). Then the triple (X°,y,(°) is a Markov process with gener¬ 
ator C. 
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2. Let (X £ ,y,( £ ) satisfy the assumptions of Theorem \2.4\ Let 

CV) = Jo ^ 9ybl(y(r))dr \l (/°° Kx,y(s))<tejb(X 0 (s), y(s))L x ° (ds, 0) (6.2) 

where the above integral is the Riemann-Stieltjes integral with respect to the increasing 
function t L x °(t , 0). Then the triple (X°, y, £°) is a Markov process with generator C. 

Proof of Theorem \2. 2\ and Theorem \2.4\ Analogous to the proof of Theorem 12.61 in the previous 

subsection. □ 

7 Remarks and further directions 

There are several interesting further questions one can ask in relation with (12.ID - (I2.2H . 

• In some sense, the original result of m and the results in this paper carry out the first 
few steps of the program contained in m, Chapter 7) when the fast process is positively 
recurrent. In particular, a result similar to Theorem 12.21 was proved by Khasminskii ([15]) 
and one sirnlar to Theorem 5 was proven again by Borodin ([5]) sharpening a result of 
Khasminskii (jl4j) that in turn was inspired by the non-rigorous result of Stratonovich 
(|28]). The next step in this program would be to study large deviation estimates of 
(X £ ,Y £ ) from their respective limit in H7j. 

• Along the same lines, it is expected that one can generalize Theorenr l2.6l to the case when 
there are several conserved quantities of that are conserved by the limiting motion in m 
Namely, if H is such a quantity and (A 0 , Y°) is the limit of (X e , Y £ ) in distribution then 
H(Y £ (t)) will converge to the constant process process H(y). uniformly on compact time 
intervals in probability. It is natural to conjecture that this process needs to be considered 
on time-scales of order e~ 2 in order to see any non-trivial behavior. The corresponding 
result for positive recurrent averaging was carried out by Borodin and Freidlin (0)- The 
result is expected to be a combination of external averaging due to the null recurrent fast 
process and internal averaging inside the over the level sets of H. In the case of Theorem 
12.61 H(y) = y and the level sets are single points, so this internal averaging does not take 
place. 

• A different direction is to replace the fast motion with a more general null recurrent 
process or a process that is only neighborhood recurrent. For example, by considering 
a situation when the fast process is driven by a Bessel process of order n € (1,2], one 
can hope to study the asymptotic behavior of a three dimensional diffusive stochastic 
dynamical system the dynamics of which is perturbed in a narrow tube. 

• The results in this paper also imply certain results on partial differential equations 
through the well known representation formulas. For example, consider the Cauchy 
problem for u £ : R 1+d x R + —> R 

d t u e = i \ip{e~ l x,y)\ 2 d xx u e + [b l {y) + b 2 {e~ l x,y)] d y u £ + 

yd d k 

+ 2 S (aa T ) i:j (e 1 x,y) d yiyj u £ + EE °ij{e 1 x,y)ip j (e 1 x,y)d xyi u £ 
i,j =1 *= 1 j =1 
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with initial condition u £ (x, y, 0) = f £ (x, y). It is well known that 


u s (x, y, t ) = E ( x ,y)fs(X £ (t), Y £ (t )) 

is the unique solution of this Cauchy problem. In this setting, if f e (x, y) = f(y) then we 
have by Lemma [3721 that u £ converges to u°(x,y,t ) = f(y(t)) uniformly on compact sets 
which solves the transport equation 

d t u° = bi(y)d y v°, u°(x , y , t) = f(y). 


On the other hand if b\ = 0, f e (x, y) = /(ex, y), then Theorem ^.dl imolies that u £ (x, y, t/e 2 ) 
converges to u°(x,y,t ) = E^^/fA 0 , Y°). Formally, this function solves 


d t v° = -a±(x,y)d xx u° + 5(x) 


/3(y)d y u° + -a(y)d yy u° 


•u°(0,x,y) =/(x,y), (7.1) 


where S(x ) is the Dirac-delta distribution at x=0, even though the solution theory of such 
an equation is non-trivial. The case when there is only a generalized drift, i.e. a = 0, 
has been studied (see e.g. [25J, but see also ®) but to our knowledge, there have been no 
successful attempts to make sense directly of (17.111 in the general case. Another approach 
is to note that Y(t) is a fractional diffusion. Indeed, in the case of e.g. a± = 1, and 
/(x, y) = /(y), we have (see [26] or the more general m 

d) /2 u° = f3(y)d y u° + ^a(y)d yy u°, u°(0, y) = f{y ), (7.2) 

where d(j is the Caputo-derivative, i.e. 

1 f°° 

e ‘ m = fa 

This latter non-local equation expresses the averaged (non-Markovian) dynamics of the 
Y -process (which is trivial when a-t = 1 and it is expected that a similar equation can be 
written down for the general case. It would also be interesting to investiate the relation 
between (EU) and m- 

Similarly, the results in this paper can be also used to study other initial-boundary 
problems through the usual representation formulas. 


Appendix A Martingale problem - SDE equivalence in the general case 

In this section we establish that there is a one to one correspondence between the solutions 
of the martingale problem associated to the operator £ and the weak solutions of the corre¬ 
sponding SDEs under very general circumstances. Even though this result is not needed for 
the specific results of this paper, we believe it is of independent interest and useful for future 
studies. For convenience, we prove the case p .|_ = p~ = 1/2 case rigorously and then outline 
the necessary changes for the case when p + and are constants but not necessarily equal. 
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Theorem A.l. Let C be the operator on Co(M 1+d ) described in Section 2 with p + = = 1/2. 

Then the distribution of a process Z(t) solves the martingale problem associated to C if and 
only if there is a Brownian motion W and a martingale V Zl such that ( Z(t ), W, V Zl ) is a weak 
solution of \2.8\) - \2.9\) . 

Proof. Using Lemma 13.51 it is straightforward to see that for every weak solution, the distri¬ 
bution of Z solves the martingale problem associated to C. 

Now we prove the opposite direction. Assume that ( X , Y) is a Martingale solution to the 
operator C given by 

1 d 

£f(x, y ) = -A(x, y)f xx (x, y) + ^ B k (x, y)f Vk (x, y) 

k =1 

whose domain consists of sufficiently regular functions subjected to the boundary conditions 

1 1 ^ \ d 

2 /s(0+, y) - ^fx( 0-, y) + ^ Pk(y)f yi ( 0, y) + - ^ a kj (y)f ykyj (y) = 0. (A.l) 


Let T>(C) be defined the same as T>(£) but without requiring its members to decay at infinity. 
This means that it consists of those continuous functions twice continuously differentiable in 
x and y with the x derivative possibly being discontinuous at x = 0. Since 

Y(t )) - f(Xm - f £f(X(s), Y(s))ds , 

J o 

is a continuous martingale for every / E T)(£), it is not hard to show that it is a local martingale 
for all / E V(C). 

First, we claim that the process Y(t) — Jq B(X(s),Y(s))ds does not move unless X(t) = 0. 
To see this, let \x\ > 5 > 0 and r = inf{f > 0 : |AT(t)| < 5}. Then by the martingale problem 
with a smooth function satisfying 


fix , y) 


yk \x\ > 5 

0 |x| < 6/2 ’ 


we have that 


Y k (r At)- 


rAt 


B k (X(s), Y(s))ds 


Jo 

is a local martingale. Applying the same procedure with shows that this local martingale has 
quadratic variation 0 and therefore does not move. The boundary conditions are not relevant 
for this stopped process because we stop the process before it reaches the boundary. To study 
the motion at the boundary, we note that for any k € {1 ,..., d}, the function 


f(x,y) = yk - \x\Pkiy) 

satisfies the boundary conditions (1A.1D . Because (X,Y) solves the martingale problem, 

Y k (t) - | X(t)\/3 k {Y(t)) - y k - \x\p k {;y) - f B k {X(s), Y ( s))ds- (A.2) 

Jo 

~ [' E ^(YisABjiXis), Y(s))ds 


l=i 
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is a local martingale. Note that because f(x, y) = x satisfies the boundary conditions, we know 
that X(t) is also a local martingale. Therefore by the Meyer-Tanaka formula, 

|X(t)| = |x| + L x (t, 0) + local martingale. 

Next, we observe that for any partition 0 = to < t\ < .. < tjj+i = t, 


\X(t)\/3 k (Y(t )) - \x\p k (y) = 

N N 

= Tj k (Y(tM\X(t i+1 )\ - |X(tj)|) + ^2 l^(^+i)l (Ph(Y(ti+ 1 )) - (3k(Y(ti ))). 
2=0 2=0 

As the partition gets finer, 


N 

lip yZ/3k(Y(ti))(\X(t i+1 )\ 

in L 2 (Q). As for the second sum, 



f3k(Y(s))L x (ds, 0) + local martingale. 


N 

^2 l^»+i)l {Pk(Y(ti+ 1)) - Pk(Y(ti))) 

i=0 

= T. |-X'(2i + i)| (y(U+i) + J B{X(s),Y(s))ds S j -/3 k (Y(ti) + J B(X(s),Y(s))ds 

where 

Y(t) = Y(t)~ [ B(Y(s))ds, 

Jo 

which does not move unless X{t) = 0 with probability one. Given any partition, P of [0,2], we 
can refine the partition (randomly) in the following manner. For any I = € P, 

1. If there exists s € / such that X(s) = 0 but X(ti + \) / 0, we split the interval, /, at the 
random time 

t = sup{s € (U,t i+ 1 ) : X(s) = 0} 

2 . Otherwise, we do not rehne I. 

Call the randomly refined partition P. This partition has the property that if [ 2 j, 2 j+ 1 ] € P, 
then either X{ti + \) = 0 or Y(fj_|_i) = Y(ti). Switching to this random family of partitions, 



N 


i+l 


lim y2 \X(t i+1 )\ ( f3 k ( Y(ti + 1) + / B(Y(s))ds] - f3 k (Y(U ) + f B(Y(s))ds 

l|P|po“ V V Jo J \ Jo 


r t d 


,d/3 k 


^|X( S )|^(Y( S )) J B i (Y( S ))ds 


/0 i=i 
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k= 1 


+ V/ fy k (X(s),Y(sm(Y{s))L x (ds, 0) + - / f xx (X(s),Y(s))d {X) s 
/,= ] 2 - y 0 

nE/ fy 1 y ] (X(s), Y(s))d (Mi , M,-) S + ^2 [ fxy(X(s ), T(s))d (X, Mp 


*,i=i i=i 

+ local martingale. 

We compare this to the martingale problem. Let / £ CpM 0 ^ 1 ) and define 


= /0,y) - *22 \x\Pk(y)f yk (0,y) - ^ M«L'(y)/yi%(°>2/) 

fc=i *j=i 

This g satisfies the boundary conditions (1A.1D and therefore 


g(X(t). Y(t)) - / Cg(X(s), Y(s))ds 

J o 

is a martingale. That is, 

d 

f(X(t),Y(t)) = f(x,y ) + ^2 (\ x (t)\Pk(Y(t))f yk (0,Y(f)) - \x\/3 k (y)f yk (0,y)) 


k =1 


d „ t 
0 


d 


J 2 L (Pk(Y(8))f yh (P, Y(s))) Bj(X(s), Y(s))ds 


k,j =1 


(A.3) 


in L 2 (fl). In this way, we have proven that 

\X(t)\p k (Y(t)) - \x\My) = r p k {Y(s))L x (ds, 0) + f V] 1*001 ^(Yis^BjiXis), Y(s))ds 

Jo Jo pt dy 3 

and by dS2D, 

M(t) :=Y(t) — y — f B{X(s),Y(s))ds - [ (3(Y(s))L x (ds,0) 

Jo Jo 

is a local martingale. By the Ito formula, for any / £ C 2 (M d+1 ), 

^ rt 

f(X(t), Y(t )) = /(x, y) + ^22 ) fy k (X(s), Y(s))B k (X(s), Y(s))ds 

1 


(A.4) 


+ ^2 (\ X(t J\ ( Yj( Y (B)fy,.yj ^ 0 , Y (t)) ~ \x\cx lJ (y)f m y, j ((). V)) 

i,j=1 

d nf O 

~ ^2 l*00|~— (aij(Y(s))fy iy .(X(s),Y(s))) B k (X(s),Y(s))ds (A.5) 

i,j,k =1 y/c 

d rt i ^ /»t 

+ ^22 fyk( X ( s ^ Y ( s V B k( x ( s ),Y(s))ds + - ^ / A(X(s),Y(s))f xx (X(s),Y(s))ds 


k=i J0 
+ local martingale. 


*d=i 


•'o 


(A.6) 
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By (IA.3I) . it follows that 


and 


\X(t)\/3 k (Y(f))fy k (Q,Y(t)) - \x\p k (y)f yk (0,y) 

d n 

- Y I Wk(Y(s))fy k ( o, TO))) Bj(X(s ), Y(s))ds 


J=1 J0 


Jo 


= I Pk(Y(s))f yk (X(s),Y(s))L x (ds,0 ) 


|X(i)|a i: ,(F(t))/ w% .(0,^(i)) - \x\aij(y)f yiyj (0, y) 


d r t 


d 


Y I \ x ( s )\g^{ a ij(Y(s))f yiyj (X(s),Y(s))) B k (X(s),Y(s))ds 


k =1 ' 
rt 


= I OLij{X{s),Y(s))fy iy AX(s),Y(s))L x (ds,0) 


ViVY 


Then by comparing (I A.4[) and (|A.5I1 we can conclude that 


(X,Mj) t = 0, 

( x )t = [ A ( x (s), Y(s))ds, 

Jo 

(Mi,Mj) t = f aij(Y(s))L x (ds, 0). 

Jo 


Finally, we define 

W(t):= f A~Y 2 (X(s),Y(s))dX(s), V X (t) := f a~ 1/2 (X(s),Y(s))dM(s). 
Jo Jo 

By Levy’s theorem, W(t ) is a one-dimensional Brownian-motion. Similarly, if 

L -1 (s) = inf{i > 0 : L x (t, 0) = s}, 


then by Levy’s theorem 

W 0 (s) := V x {L~\s)) 

is a d-dimensional Wiener process, and V x (t ) = Wo{L x (t , 0)). Finally, we can conclude that 
(A(f), Y(t j) is the weak solution of 

f dX{t) = V A(X(t), Y(t))dW(t), fA ^ 

| dY(t) = B(X(t), Y(t))dt + /3(Y(t))L x (dt, 0) + \/a(Y{t))dV x {t). 

□ 
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Remark A.2. In the case when p+(y) and p~(,y) are not necessarily equal, one direction is 
still given by Lemma 1 3 . 51 The main difference in proving the other direction is that now 
X(t ) is not a local martingale. Rather, by applying the martingale problem to the function 
fix, y) = x — (p+(y) — P-{y))\x\ (which requires sufficient regularity from p±), we see that 

X(t) - ( p+(Y(t )) - p-(Y(t))\X(t)\ (A. 8 ) 

is a local martingale. The difficulty lies in the fact that we do not a priori know that X(t) 
is a semimartingale and therefore we cannot immediately apply the Tanaka-Meyer formula. 
However, let Tq = 0, Og = inf{t > 0|A(t) = 0} and recursively define 

rf = inf{t > crf_ 1 ||A(t)| = <5}, of = inf{£ > r/| A(f) = 0} 

and we can write 


X(t) - A(0) = J2 [X<t A af) 

i =0 


X(t A 4) 


+ E[^(* Ar/) 

2=1 


X(t A 4.-\) 


It can be shown similarly as in Lemma 2.2 in m that as 5 4- 0, the first sum converges to a 
local martingale while the second sum converges to a monotone process that only changes when 
X(i) = 0. More precisely, this second term converges to f^(p + (Y(s)) — p-(Y(s)))dL x (s). This 
shows that X is a semimartingale and II A. 80 implies that X(t) — fg(p + (Y(s)) — p-(Y(s)))dL x (s) 
is a local martingale. The rest of the proof goes through with minor modifications. 
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